Quantum fluctuations, pseudogap, and the T = superfluid density in strongly 

correlated d-wave superconductors 
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I study the effect of Coulomb interaction on superconducting order in a d-wave lattice super- 
conductor at T — by considering the superconducting saddle point in the two dimensional t-J-U 
model with a repulsion U. The theory of low-energy superconducting phase fluctuations around 
this saddle point is derived in terms of the effective hard-core bosons (representing the density of 
spin- up electrons and the phase of the order parameter), interacting with the fluctuating density 
of spin-down electrons. Whereas the saddle-point value of the superconducting gap is found to 
continuously increase towards half filling, the phase stiffness at T = has a maximum, and then 
decreases with further underdoping. Right at half filling the phase stiffness vanishes for large U. 
This argues that the pseudogap phenomenon of the type observed in cuprates is in principle possi- 
ble without a development of any competing order, purely as a result of growing correlations in the 
superconducting state. Implications for the finite temperature superconducting transition and the 
effects of static disorder are discussed qualitatively. 



I. INTRODUCTION 

Pseudogap phenomenon has become one of the hall- 
marks of high temperature superconductivity [lj : while 
the suppression of the single particle density of states is 
observed in many quantities at a high temperature T*, 
which appears to increase towards half filling, the super- 
conducting transition temperature T c , together with the 
T = superfluid density, at the same time continuously 
approaches zero Q. Being in dramatic contrast to the 
well understood behavior of the standard BCS supercon- 
ductors 3] , this puzzling phenomenon has prompted dif- 
ferent explanations. These may be grouped into at least 
two conceptually separate camps. The first group postu- 
lates development of a second order parameter in the un- 
derdoped region, which competes with superconductivity 
and suppresses its transition temperature 0- Q> 

1 1 [HI ■ The second interprets T* as a crossover tem- 
perature at which the Cooper pairs, loosely speaking, be- 
gin to form, and the lower T c as the point where the phase 
coherence finally sets in [Tl|. 0, [3. The suppression 
of T c is then attributed to the gradual loss of carriers 
near the insulating state at half filling. There is a signif- 
icant amount of possible empirical support for the latter 
point of view: d-wave symmetry of the pseudo gap |14| . 
specific heat measurements 0] , heat transport Jlfj , mi- 
crowave conductivity 0] , and the Nernst effect |lSf may 
all be understood as directly or indirectly supporting the 
proposed superconducting origin of the pseudogap tem- 
perature. Some of the same measurements, however, may 
also be understood within a competing theory, and the 
physics of underdoped high temperature superconductors 
at present still remains controversial. 

In this paper it is argued that the pseudogap phe- 
nomenon may in principle arise purely from strong 
Coulomb interactions in a d-wave superconductor and 
sufficiently near half filling, without any competing or- 
der. While a competing order is possible, and indeed in 
some materials may even be likely in the underdoped re- 



gion, it seems not to be required simply by the existence 
of pseudogap. Furthermore, this is shown for a supercon- 
ductor with only a weak attraction in the d-wave channel, 
so the mechanism behind the pseudogap considered here 
is different than in the real-space pai ring approaches ex- 
tensively discussed in the literature [Fj. Deriving 
from Coulomb repulsion, it is similar in spirit, although 
still different in detail, to the one in the RVB theory [2(J, 

m, p. 

We consider what is probably the simplest model which 
contains the relevant physics, the t-J-U model in Eq. (1), 
with the standard exchange term rewritten as the pair- 
ing interaction and with a strong repulsion, and study 
the evolution of the T = superconducting gap and the 
phase stiffness with doping. Independently of its mag- 
nitude, the repulsion does not affect the dependence of 
the mean-field d-wave superconducting (dSC) gap on the 
number of particles, as it only shifts the (unphysical) bare 
value of the chemical potential. As a result, the saddle 
point value of the superconducting gap (Figure 1) con- 
tinuously increases towards half filling, as found in other 
similar calculations |23|. It is intuitively clear however, 
that for a large repulsion U(x) the superconducting order 
should eventually be weakened near half filling, since the 
Cooper pairs, although formed, will have little space left 
to move. It is less obvious in which exact fashion should 
this occur; for example, should the suppression be grad- 
ual or abrupt 1241. or what its effect on quasiparticles 
should be (23, |26|. It is possible to make the physics 
behind this 'jamming' effect quite explicit by introducing 
certain collective coordinates to describe the (quantum) 
phase fluctuations around the dSC saddle point: two den- 
sities of spin-up and spin-down electrons, and the phase 
of the superconducting order parameter. The crucial fea- 
ture of our representation is that only one of the electron 
densities (say of spin-up) becomes the conjugate variable 
to the superconducting phase. This enables one to un- 
derstand quantum (T — 0) phase fluctuations around the 
dSC saddle point within the theory of interacting bosons 
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(representing jointly the density of spin-up electrons and 
the order parameter phase) moving in the fluctuating and 
interacting background provided by the density of the 
remaining spin-down electrons. A self-consistent calcu- 
lation of the T — phase stiffness in such an effective 
bosonic system yields then generically a non-monotonic 
behavior with doping (Figure 2). For U large enough 
the phase stiffness at half filling vanishes, while the su- 
perconducting gap remains perfectly finite, and in fact 
large. 

The superconducting T c is then non-monotonic be- 
cause it is determined by the smaller of the two character- 
istic temperatures [l^ . one set by the superconducting 
gap, and the other by T = superfluid density. It is 
not the bare superfluid density that vanishes near half- 
filling, however, but only the physical (renormalized) one 
that becomes reduced by quantum fluctuations. The bare 
one, indeed, only increases towards half filling, since in a 
weakly coupled superconductor it is proportional to the 
electron density. Furthermore, in contrast to T c T = 
superfluid density should continue to grow as doping is in- 
creased through the optimal doping, as observed in num- 
ber of experiments |2jfl . The characteristic temperature 
dependence of the superfluid density under the condi- 
tions of a large d- wave gap and a low T c is discussed in 
a separate publication |28j . 

In the previous work on the d-wave quasiparticles [2(J 
it was assumed that underdoping, known from experience 
to be detrimental to superfluid density, could be viewed 
as inducing quantum phase fluctuations into an other- 
wise standard d-wave BCS-like superconducting ground 
state. Present results support that hypothesis. It was 
also shown that such quantum phase fluctuations lead 
generically to spin response which agrees in many de- 
tails with the observations in underdoped Y-Ba-Cu-0 
(YBCO) |29j. Together with the present results, this 
lends additional credence to the idea that the underdoped 
cuprates may be understood as strongly (quantum) phase 
fluctuating d-wave superconductors. The effect of low 
phase stiffness on charge degrees of freedom will be dis- 
cussed in a future publication |30j| . 

The paper is organized as follows. In sec. II the 
t — J — U Hamiltonian is introduced, and the supercon- 
ducting saddle point is discussed in the functional inte- 
gral formulation. The T — theory for the phase fluctu- 
ations around the saddle point is derived in sec. III. The 
effective bosonic theory is presented in sec. IV, together 
with the self-consistent calculation of the T — super- 
fluid density. The discussion of the results in relation to 
the physics of cuprates is given in sec. V. 



II. MODEL AND THE SUPERCONDUCTING 
SADDLE POINT 

We begin by defining the t-J-U Hamiltonian. Simi- 
lar models have been often invoked in the past to study 
strongly correlated electrons [U, [S^, [33. Consider 
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FIG. 1: The saddle point of the gap Ao/t vs. doping x — 1—n, 
from Eqs. (4) and (5). 



the imaginary time quantum mechanical action S = 
J o drL(r), P = 1/T, and L(r) = E x ,^U^r)(d T - 
/j,)c <t (x,t) + H(t), with c a (x,T) being the standard 
Grassman electronic variables, and the Hamiltonian 

H{r) = -t ]T 4(x,t)c ct (x',t)+(I) 

{x,x') ,(T — ± 

^(n + (x,r) +n^{x,T)) U ^ X ^ — {n + {x',r) +u_(x' 1 t)) 

X,X f 

+J ^2 [S(x,t) ■ S(x',t) - -^n CT (a;,T)n CT '(x',T)]. 

{x.x') cr.a' 

Here, S(x,t) = (1/2) J2 a ,p c ai x i T)a a pcp{x 7 t) and 
n (J {x 1 t) = c£(x, t)c ct (x, t) are the standard electron spin 
and particle densities, U(x), J > 0, and x labels sites of a 
two dimensional quadratic lattice. When U(x) = U5 XtX i 
with U — ► 00 the model becomes equivalent to the stan- 
dard t- J model as derivable from the underlying Hubbard 
model, since double occupancy then becomes completely 
suppressed. Here we will relax this condition, and con- 
sider a general, and finite repulsive interaction U(x). We 
will however, be particularly interested in the regime of 
parameters U(x) 3> t > J, which corresponds to a weakly 
coupled, but a strongly correlated superconductor. 

It is presently controversial [3^] , |3^] whether the 
canonical t-J model indeed supports the d-wave super- 
conducting ground state. Although interesting, this issue 
will be of little concern to us here. Our philosophy 36], 
[25|. ptl shared by number of recent works |22j . |37j. is 
to start from the experimental fact that cuprates are d- 
wave superconductors at least for a range of dopings [38| , 
and take that as the point of departure for further explo- 
ration of the phase diagram. We will therefore regard the 
Hamiltonian (I) merely as a useful tool for introducing 
correlations into the postulated superconducting ground 
state. 

The last (~ J) term in the Hamiltonian can be rewrit- 
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ten to explicate its pairing nature as 



J 



Hj = -- B^x,x',T')B(x,x',r), 



(2) 



with B(x,x') — C-(x' ,t)c+{x,t) — c+(x' \t)c-(x,t) an- 
nihilating a singlet on a pair of neighboring sites. In the 
momentum space, Hj is a sum of the d-wave and the 
extended s-wave pairing terms [23 ]. In what follows we 
consider a purely d-wave saddle point, and neglect com- 
pletely a possible s-wave component. Being small at the 
Fermi surface near half filing, the amplitude of the s-wave 
order parameter is expected to be completely suppressed 
by the dominate d-wave component 39]. 

We start by decoupling the Hj term using a complex 
Hubbard- Stratonovich field A(x,x' ,t), and by introduc- 
ing two density variables p (7 (x, r), u = ± as 

L(t) = 2J c ff( x > r )[(^T - V + ijo-( x )) s x,x' -t]c„(x' ,t){2>) 

(x,X f ) ,(T 



-i^23v(x,T)p a (x,T) 

x.a 



2 p ( x > T ) U ( X ~ x ')p( x> , T ) 



+ £ [^\A(x,x',T)\ 2 + A(x,x',r)B^ a 



H.c 



[x,x') 



The variables ja(x,r) are the Lagrange multipliers that 
enforce the constraints p a {x,T) — n a (x,T) 40], and 
p(x, t) = p+(x,r) + p- (x, t). The d-wave superconduct- 
ing saddle point is as usual given by A(x, x', r) = +A if 
x' = x + 1 and A(x, x',t) = — Aq if x' = x + 2, where 1 
and 2 are the unit lattice vectors, and with the amplitude 
Ao 7^ at T = being determined by the standard BCS 
gap equation 



1 = J 



d 2 k 



<t>m 



(27T) 



(e(fc)-A*) a + |Ao| a ^(fc) 



(4) 



together with the number equation 



d 2 k 



e(k) — p 



(e(fc)-/i) 2 + |Ao| 2 ^) 



(5) 



n represents the number of electrons per site, related to 
doping x as n = 1 — x. Here e(k) = —2t(cos(k x ) + 
cos(fcj,)), (j)d(k) = cos(fc x ) — cos(fcj,), and p — p — nll(q = 
0). We assumed a non-magnetic saddle point at which 
p+ = P- = n/2, and j+ = j- = —inU(q — 0) |4lj . 
Our goal here is to consider quantum fluctuations result- 
ing purely from the electron repulsion, i. e. without any 
explicit competing order. This is not to say that some 
other saddle point may not become stable as the elec- 
tron density or the interaction U is varied; indeed, in the 
particular t-J-U model this would be expected for large 
U . Our concern here, however, is not the phase diagram 
of the specific t-J-U model [33 , but the broader issue of 



the evolution of the postulated superconducting ground 
state in presence of strong repulsion. 

An example of a numerical solution of the saddle-point 
Eqs. (4)-(5) is presented in Figure 1, for parameters 
t = 150meV and J = lOQmeV. Even with such a rel- 
atively large J the superconductor is effectively still in 
the weak-coupling regime, as evidenced by the chemical 
potential, for example, changing very little by the open- 
ing of the superconducting gap. In the rest of the paper 
we will therefore be concerned exclusively with such a 
weakly-coupled superconducting state. Figure 1 should 
of course be understood only as an illustration, as a more 
realistic calculation should include at least a sizable next- 
nearest- neighbor hopping t' . It does demonstrate, how- 
ever, two expected and generic features: a) the gap at a 
given electron density is unaffected by the repulsion U (x), 
which only shifts the bare value of the chemical potential 
p, b) the saddle-point value of the gap is a monoton- 
ically increasing function of the electron density. The 
mean-field analysis by itself would therefore suggest an 
increasing superconducting T c towards half filling, in bla- 
tant contradiction to the generic behavior of underdoped 
cuprates. This conclusion will be overturned however by 
the inclusion of quantum phase fluctuations introduced 
by the large repulsion, as discussed next. 



III. QUANTUM FLUCTUATIONS 

To this end, consider the fluctuations of the phase of 
the order parameter, A(x,x + a, r) = (— ) a Aoe~ ie ( x ' T \ 
where a = 1, 2. To reduce the algebraic complexity we 
will neglect the fluctuations of the amplitude of the or- 
der parameter, which should be justified at T — 0. We 
are also assuming a single phase 9(x, r) for both links 
(x, x + i) and (x,x + j) emanating from the site x. In 
doing so we are neglecting the fluctuations towards the 
possible s-wave component of the superconducting order, 
in accord with the expectation that it is suppressed by 
the dominate d-wave order [39| . 

As usual, at this stage one would like to integrate out 
fermions, to be left with an effective action for the col- 
lective variables. To perform this step, we first introduce 
new Grassman variables as c\_(x,t) = e iei - x,T ' l a ! + (x, r), 

c + (x,t) = e~ l0( - x ' T ^a + (x, t), c_(x,t) = o,_(x,t), and 
C-(x, t) = a-(x,r). Note that in contrast to previous 
works [42]] , [43] ] , the entire superconducting phase here is 
'absorbed' into particles of a single spin-projection, ar- 
bitrarily chosen to be up [44|. This guarantees that for 
all possible configurations of the phase, including vor- 
tices, new Grassman variables a and a' are single valued 
and satisfy the standard fermionic boundary conditions 
at t = j3. Although the above change of variables is not 
unique in accomplishing this it probably is the simplest. 

Next, as usual, introduce the deviations from the 
saddle-point values of the fields as ij a {x, t) — Sj^^x, r) + 
nU(q = 0), p a {x, t) = 5p a (x, t) + n/2, a — ±, and then 
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shift Sj + (x, t) — id T 9(x, r) — > 5j+(x, r). The Lagrangian 
then becomes L(r) = L a (r) + LJt), with 

£a( T ) = X! t( a; ) T )[( a r _ / i + ^'er( a! J T ))^>i=' ( 6 ) 
(x,x f ) ,<J 

£ ±A [ot(ar, r)al (x', r) - e -*W^)-*(*V)) 

(x,x') 

a'L(x,T)a t + (x',t) + H.c], 
L p (t) = -^25j a (x,T)p a (x,T) + (7) 

x,a 

\ ^2 Sp ( x > T ^ U ( x ~ v) 5 P(v> T )~ i ^2p+ i x > r)d T 6{x, r). 

x.y x 

The shift of one of the Lagrange multipliers serves to iso- 
late the imaginary-time derivative of the phase into the 
last term in L p , and thus promote the density of spin- up 
electrons p+ to the status of conjugate variable to the 
superconducting phase. It is important to note that it 
is not (p + (x,t) + p-(x,r))/2 that plays the role of the 
conjugate variable, as it would have been obtained if the 
two projections of spin were treated equally, and half of 
the superconducting phase absorbed into each species. 
That, often invoked transformation is allowed only for 
the trivial boundary condition Q(x,($) = 9(x,0) + 2irn, 
with n = 0, and with the vortices in the phase being for- 
bidden. For this trivial boundary condition the last term 
in Eq. 7 becomes independent of the average particle 
density, and any effect of the proximity to commensu- 
rate density, which is our main subject, is lost. To have 
any commensuration effects on the superfiuid response 
it is thus paramount to allow for the non-trivial (in the 
above sense) phase configurations. 

In principle the fermions can now be integrated out. 
The resulting action will be a functional of the phase 
9{x, t) and the two densities p±(x, r) of spin-up and spin- 
down electrons. The fact that p+ and 9 form a pair of 
conjugate variables suggests a particular change of vari- 
ables, so that the action becomes a functional of the 
bosonic variable \P = y/p+e l9 , \&*, and This is the 
central idea of this paper. In the rest we will try to ap- 
proximately determine the functional S[^, \&*, and 
then use it to compute the stiffness for the superconduct- 
ing phase fluctuations at T = 0. Alternatively, the lattice 
action in Eqs. 6 and 7 can be studied numerically. 

In practice however, the integration over fermions is 
performable analytically only for small phase gradients 
and a small 5j a . After a straightforward but a rather 
involved algebra one finds S = Sf + Jq ' dTL p (r), with 
the contribution from the fermion determinant 

sp= \S S^ 1 ® v)m v) + (8) 

\ f dr^2Sj a (x,r) + 0(R% 



where Fft(q, v) = (Sj + (q, v), —i9(q, v), 5j-(q, v)), and the 
3x3 matrix M is symmetric with the elements 

/i 2 
— ^G(k,u)G(k+q,w+v), 

(9) 

M 13 (q,v) = -J ^^F(k,u)F(k + q,u; + v), (10) 

M 12 ( 9 » = J^-G{k + q,u> + v)[{e{k + q) (11) 

-e{k))G(k, u) + (A(k + q)- A(k))F(k, uj)} 
+ {A(k + q)- A(k))F(k + q,UJ + u)G{k, u), 

M 23 (q,u) = J *^F$+?,w + v)[-(e(k + $ (12) 

-e{k))F[k, uj) + (A(k +q)- A(k))G{-k, -«)] 
+ (A(k + q)- A(k))G(-k -q,-cj- v)F{k, w), 

and finally 

M 22 (g» =(13) 

/ [e{% +q) ~ e &w G & + + ") 

+2[e(k + q)- e(k)](A(k + q) - A(fc))[G(fc + q, uj + v) 
F(k, u) + G(k, uj)F(k + q, uj + v)\ + 2[A{k + q) - A{k)} 2 
[F(k, uj)F{k + q,uj + is)- G(-k, -uj)G{k + q,u + v)] 

+2[e{k + q)-e{k)]G{k,w) 
-2[A(-k + q)- A(-k)]F(k,uj), 

withG(k,uj) = (-iuj~e(k)+p)/[uj 2 +A 2 (k) + (e(k)-p) 2 }, 
and F(k,u) = A{k,uj)/[uj 2 + A 2 (k) + (e(fc) - p) 2 }, and 
A(fc) = A (j) d (k). Note that Eqs. (9)-(13) are similar, 
but not identical to those that can be found in (4^, for 
example. This is due to a different change of fermionic 
variables employed there which treats spin-up and spin- 
down symmetrically, as discussed earlier. In particular, 
we have M12 ^ M23, and as a result the two densities 
are in principle coupled differently to the gradients of 
the phase. This serves to compensate for the asymmetry 
between p + and p_ in L p , and to insure that if solved 
exactly, the theory would lead to equal average densities 
of up and down particles. 

The elements of M(q, v) can be easily evaluated 
for a weakly coupled superconductor with Ao <C t. 
Assuming a spherical Fermi surface for convenience, 
we find M n (q,v) = -JV/4 + 0{q 2 ,v 2 ), M 13 (q,v) = 
-N/4 + 0(q 2 ,v 2 ), M 12 (q,is) = M 23 (g» = 0(q 2 ), and 
M 22 (q, v) = (e F /47r)(l + 0((A Q /e F ) 2 ))q 2 + 0{q\ q 2 v 2 ), 
where M is the density of states at the Fermi level e F 
(M = m* /tt for the spherical dispersion e(fc) = k 2 /2m*, 
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with m* — l/(2ta 2 ) being an effective eiectron mass, and 
a the lattice spacing). 

At long distances (q, v — > 0) one can therefore neglect 
the matrix elements Mi 2 and M23 that couple the fluctu- 
ating phase to the Lagrange multipliers. The integration 
over Sj± then gives the action in terms of the densities 
and the phase to be 



S = 



d 2 qdv ep 2 



(2tt) 3 l 8tt 



{^e(q,v)9(-q,-p)+ (14) 



■(dp + (g,i/) + op_(g,i/)) 



2(M 2 3 (g>)-M 1 2 1 («?->)) 



M 13 (q,v) 



M* 3 (q,v)- M^(q,v) 



Sp+(?,u)Sp-(-q, -v) + 



-U{q)6p(q, v)5p(-q, -v) + v6(q, i/)p + (-q, -v)}, 

where we retained only the leading terms in the expan- 
sion of S in powers of phases, densities, and their deriva- 
tives. The bare phase stiffness ~ e.p = 2imt for weak 
coupling, and in general will be proportional to the total 
kinetic energy of electrons in the superconducting state. 

The integration over fermions has introduced addi- 
tional quasi-interactions between the electron densities 
in Eq. 14, on top of the repulsion U{q): 



M u (q,is) 



~ (Mi 2 3 (?, v) - Mi 2 i(?, v)) 
1 

2(Mi 3 (g»-Mn(q»)' 



(15) 



and 



V-+(q,u) 



M 13 (q,v) 



2{Mi 3 {q,v) - M u %v)) 



(16) 



The quasi-interaction is repulsive between the same and 
attractive between the opposite spin densities. The ap- 
proximate expressions in Eqs. (15) and (16) are valid at 
small q. Here we utilized the fact that (for both d-wave 
and s-wave superconducting state) at T = 



Mi3(0,i/) = Mu(0,i/), 



(17) 



exactly. This is easily established using the definitions of 
the matrix elements and performing the frequency inte- 
gral. Physically, this simply expresses the fact that the 
uniform spin susceptibility in a superconductor with sin- 
glet pairing vanishes at T — 0. This is probably most 
easily seen by 'rotating' our density variables into total 
density and magnetization, as p = p++p-, m = p + — p_, 
and then reading off the [random-phase approximation 
(RPA)] spin and charge susceptibility in the supercon- 
ductor from Eq. (14): 



X s 1 (<?» 



1 



2(M 13 (g»-M 1 i(g»)' 



(18) 



X - 1 (q,v) = U(q) 



1 



2(Mi 3 (g»+Mu(g»)' 



(19) 



Eq. (17) therefore implies a vanishing uniform spin sus- 
ceptibility and a finite compressibility (for short-range 
interaction) in a singlet superconductor at T = 0. Note 
also that for non-interacting electrons (U = J = 0) one 
recovers \s = Xc — -2Mu(g, v), which is then the fa- 
miliar Lindhard response. 

We may gain further insight into the nature of quasi- 
interactions V++ and V-\ by computing them at small q. 

In the first approximation one may neglect the anisotropy 
of the d-wave gap and perform the calculation for the 
simpler s-wave superconductor. The result is 



M 13 (q,v) - M n (q,v) 



67rA§ 



q 2 (l + 0(v 2 )). (20) 



We checked numerically that the result for d-wave su- 
perconductor is essentially the same, except for the 
anisotropy in the coefficient of the g 2 -term. Taking the 
d-wave nature of the gap into account yields a maximum 
of Mi 3 — M11 at a finite q that connects the gap nodes, 
so that the Stoner instability of a metal in a dSC be- 
comes replaced by an analogous RPA instability towards 
the spin density wave 0- For the present purposes it 
suffices to note that the quasi-interactions induced by 
the fermion integration are long-ranged, but weak, for a 
weakly coupled superconductor, i. e. V aa < ~ Ag/e^. 
In particular, this means that the electrons with the op- 
posite spin, although paired in the momentum space by 
construction, may to a good approximation be considered 
as independent in real space, and hence to be interacting 
primarily via U(x). 



IV. EFFECTIVE BOSONIC HAMILTONIAN 
AND THE PHASE STIFFNESS 

If one would integrate out the Gaussian density fluc- 
tuations in Eq. 14 the remaining action for the phase, 
apart from the linear time derivative term, describes the 
standard collective mode in a weakly coupled supercon- 
ductor. To capture the effects of nearly commensurate 
density on the superfluid response one needs to go be- 
yond this hydrodynamic regime What is needed 
is the bosonic lattice action S[$>, \&*, p_], which would 
correspond to Eqs. 6 and 7, and which at long length 
scales would reduce to Eq. 14. The simplest candidate 
satisfying these requirements is: 



iB(r)=^**(a;,r)S T *(ar,r)- 

X 

t B J2 ^*(x,t)^(x',t) 

(x,x') 

■\ J dr'^2mx,r)\ 2 (V ++ (x-x',r-r') 



(21) 



+U(x - x')5(t - t'))\^(x' ,t')\ 2 
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place the realistic Coulomb interaction by the hard-core 
lepuision DCLwecn me paiucies. -Lo niaKC ine calcula- 
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tion a little more interesting however, we will assume the 
hard-core repulsion only between the same spin particles 
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[49j] , and allow for some finite on-site repulsion JJ between 
the particles with opposite spin. The single parameter U 


. 4 




then effectively measures the degree of correlations at the 
superconducting saddle point. For U — oo then all the 


0.2 




particles interact via hard-core repulsion. 







The phase stiffness from Eq. (21) is therefore 



0.2 



0.4 



. 6 



0.8 



FIG. 2: The phase stiffness K/Ats in the effective bosonic 
theory for phase fluctuations for (7 = 0, U/4tt = 1.8, U /Ats = 
5, and tsjU = (top to bottom) vs. doping x = 1 — n. 



+ J <1t'^\^{x.t)\ 2 (U{x - x')5{t - t 1 ) 

x,x' 

+V + _(x-x',t-t'))p-(x',t') 

X 

\ J dr> E p-( x i t )( v --( x - x> > t - t ') 

ie, a; 7 

+U(x-x')6(t-t')) P -(x',t'), 



with *(x,r) = ^p + (x,T)e t0< - x > T \ t B = e F /(4nir) and 
p,B = xJ 1 (0,0)n/2. The reader is invited to check that 
when expanded in powers of the density and phase and 
their derivatives [47] . the leading order terms in L B in- 
deed reproduce the Eq. (14). 

The action in Eq. (21) represents bosonic particles 
hopping on a lattice and interacting with a fluctuating 
background provided by the density of spin-down elec- 
trons. Although the amplitude of our bosonic field is 
proportional to the density of spin-up electrons only, its 
phase is still the full superconducting phase and there- 
fore the bosons have the electromagnetic charge 2e. Since 
ep ~ tn, hopping amplitude for the bosons is is = 
eF/(4n7r) ~ t, and is approximately doping independent. 
Since the change of fermionic variables leading to Eq. (7) 
broke the spin up-down symmetry, although the chemi- 
cal potential fi B is the same for both the bosons and the 
spin-down particles in Lb, two densities in our approx- 
imation will not automatically be equal as they should. 
One may easily correct this by allowing for two different 
chemical potentials p + and and adjust them so that 
the average densities of bosons and spin-down electrons 
become equal, as of course they are in the full theory 
in Eqs. 6 and 7. Furthermore, assuming the Coulomb 
interaction U(x) ~ 1/x, for a weakly coupled supercon- 
ductor we may neglect the quasi-interactions V-\ and 

V ++ , which are weak, and also less singular than U(x) 
at x = |4^. To simplify the calculation that follows, 
after neglecting the quasi-interactions we may further re- 



K = t B \{9(x,T))\* 



(22) 



Within the local approximation discussed above we 
may compute K in the effective bosonic action in 
Eq. (21) by decoupling the hopping term with a 
Hubbard- Stratonovich field $(x, t), where ($(x, t)) = 
2t B ('i'(x, t)), and by making yet another saddle-point 
approximation in the resulting action. Assuming that at 
such a saddle point both $(x, r) and p_(x,r) are inde- 
pendent of the imaginary time, the 'mean-field' bosonic 
free energy becomes (see Appendix) 

F b ,mf = 2E**( a 0(*5 1 )x,«'*0c / ) (23) 



x,x' 



1 rl 



-T^lnTrExp - — [-®(x)&(x) + -$*(x)#(x) 

X 

-H+&(x)i>(x) - /x_p_(x) + (x)i>(x)p-(x)], 

where the trace is to be evaluated over four possible 
states: empty site |0), boson ^(x^O), spin-down elec- 
tron |— ), and boson and spin-down electron ^(x)!— ). 
Since we assumed a static saddle point, the temperature 
T appears explicitly in F B , and at the end we need to 
impose the limit T — ► 0. $(x) takes the value that min- 
imizes Fb.mf, an d the chemical potentials fi± are to be 
tuned so that the densities of bosons and of spin down 
electrons are both equal to half the electron density: 



dF, 



B,MF 



<9$(x) 



= 0, 



n _ 1 dF BjM F 
2 ~ M dp+ ' 

n 1 d Fb.mf 

2 ~~ M Mdp- ' 

Evaluating the trace in F b ^mf yields: 

^<Z>*(x)(tz% tX ,$(x') 

x.x' 

4 

Mp+ 



F, 



B .MF 



(24) 



(25) 



(26) 



(27) 
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with 



DISCUSSION 



A 1)2 = ±V/4 + l$(*)l 2 /2T, 



(28) 



and 



A 3 , 4 = [2/i- - C/ ± ^(U-^ + Y + \^{x)\^]/2T. (29) 

M is the number of lattice sites. At T = one can 
therefore neglect the terms containing A2 and A4, to find 

_L_ _ ((^ + l^l 2 )- 1 / 2 ^ 1 + ((U - M +) 2 + \^\ 2 )- 1 ' 2 e^ 

(30) 

(31) 
(32) 



n _ 1, y^+lfl 2 y(t/-M+) 2 +|t>l 2 
2 2 l e A i+e A 3 



n 
2 



e Ai _|_ e A 3 



where we have also assumed a uniform $(2;) = $ 7^ 0. 

The mean-field Eqs. (30)-(32) can be solved for any U, 
but we focus here only on the limits of large and small 
interaction. When [7 = 0, one finds 

K= — =n{2-n)4t B + 0{U). (33) 

The phase stiffness of the non-interacting (U = 0) weakly 
coupled dSC is a monotonically decreasing function of 
doping x = 1 — n, behaving in essentially the same way 
as the saddle-point amplitude Ao, except for being much 
larger. This is the familiar non-correlated BCS limit. In 
the strongly interacting case lifts 3> 1, on the other 
hand, 



— = 2n(l - n) + 2n 2 (l - n) — 



(34) 



As n — > the superfluid density, of course, vanishes for 
any U . For strong interaction, however, as n — » 1 , 
if — > for U/4t B > 2. The order parameter is a de- 
creasing function of U, and its evolution with interaction 
is presented on Fig. 2. 

Physical mechanism behind the suppression of the 
phase stiffness near half filling is simple: although the 
number of bosons per site near n = 1 is only ~ 1/2, 
the remaining sites are largely occupied with spin-down 
particles. For U/Atb > 2 and right at half filling this 
completely blocks the motion of bosons and brings the 
stiffness to zero. For a small U, however, it is advan- 
tageous for bosons to be in the superfluid state even at 
half-filling, since there is only a small energy cost to hop 
through a site occupied with a spin down electron. At 
U = 0, of course, spin-down electrons are (almost) invis- 
ible to bosons (with neglected), and the stiffness at 

half filling is actually the largest. 



Our main result is that for a strong repulsion U, the 
T = superfluid density becomes strongly suppressed 
near half filling, and in fact vanishes right at for large 
U . This is demonstrated for a weakly-coupled supercon- 
ductor, in which the electrons constituting a Cooper pair 
may be considered essentially uncorrelated in real space. 
The gist of the method is the introduction of two sepa- 
rate density fields for spin up and down electrons, only 
one of which is grouped together with the superconduct- 
ing phase into a bosonic variable, while the other pro- 
vides the interacting fluctuating background. Vanishing 
of the phase stiffness towards half filling is then obtained 
without invoking any competing order parameters or the 
real-space pairing. 

That the found behavior of the superfluid density im- 
plies a pseudogap behavior and the nonmonotonic super- 
conducting T c with doping is qualitatively seen as follows. 
In principle, the superconducting T c is determined by the 
lower of two characteristic temperatures, the mean-field 
Tbcs oc Ao, and Tg oc K at which the phase coherence 
would dissapear even with the amplitude Ao held fixed. 
Since the characteristic energy scale for Tg is the hopping 
parameter t, and for a weakly coupled superconductor 
Ao <C t, for a weak interaction U Tg 3> TbcSi an d it is 
the Tbcs that determines T c For strong U however, 
this will remain true only above some doping x op t- For 
the parameters used in Fig. 1 x opt < 0.5, since atiw 0.5 
the gap becomes rather small while the superfluid den- 
sity there peaks. For x -C x op t Tg -C Tbcs and it is Tg 
that will determine T c . One can still see the remnant of 
Tbcs a $ the crossover temperature where the density of 
states becomes suppressed, just as it would be in a true 
superconductor. This divorce of Tbcs and the true su- 
perconducting T c is an unavoidable consequence of strong 
repulsion, and may be taken as an operating definition 
of a strongly correlated d-wave superconductor. 

Our T = calculation has another interesting conse- 
quence worth mentioning. While a strong U implies the 
existence of the optimal doping x op t at which T c peaks, 
at T = there is nothing particular happening at that 
doping at all. In fact, as doping is increased above Xopt, 
T = superfluid density should continue to increase. 
This is precisely what is generally observed in high tem- 
perature superconductors |27| . Experimentally, the in- 
crease of superfluid density is continued until a higher 
doping x' opt , above which the superfluid density again 
starts decreasing. We believe this effect is likely to be 
due to disorder, about which few words therefore need to 
be said next. 

We found that without any disorder, and with all com- 
peting orders excluded, T = superfluid density in a 
strongly correlated dSC is suppressed near half filling, 
but always remains finite at finite doping. Experimen- 
tally on the other hand, it appears that the superfluidity 
in high temperature superconductors is already lost be- 
low some critical doping x u « 0.05. The addition of dis- 
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order to our bosonic description of quantum phase fluctu- 
ations in underdoped regime would be expected to lead to 
precisely such a behavior, as suggested by the studies of 
superfluid-insulator transition |50|. Recalling that elec- 
trons, and consequently our effective bosons, in reality 
interact via long-range Coulomb interaction suggests the 
universality class of two-dimensional dirty bosons with 
z = 1 , and v « 1 |5f| , in accord with a number of differ- 
ent experiments in underdoped cuprates [H^ ■ In the over- 
doped regime, on the other hand, A becomes small and 
the pair breaking effect of disorder will bring it, together 
with T c and the superfluid density, to zero at another 
critical doping x |53j. The precise mechanism of this 
suppression seems also in agreement with the present, 
albeit somewhat limited, experimental results on over- 
doped cuprates |27|, This would also explain why T = 
superfluid density starts to decrease for x > x' t , where 
it would still be increasing according to our picture, if 
there was no disorder. 

Although in the present work we were concerned with 
T = 0, our picture of the quantum superconductor- 
insulator transition suggests that the universality class 
of the transition at T = T c in underdoped regime will 
be of Berezinskii-Kosterlitz-Thouless (BKT) type. In- 
deed, under the assumption of a finite superconducting 
gap, vortex unbinding becomes the only known mecha- 
nism for the loss of phase coherence in two dimensions. 
Possible BKT nature of the transition is supported by 
the observed large Nernst effect above T c in underdoped 
cuprates [Tsj . and the measurements of the microwave 
conductivity ^tJ- K was argued recently however, 
that having the BKT transition well below T* set by a 
large amplitude Ao requires a competing order param- 
eter. The argument is as follows: the condensation en- 
ergy per unit area of the BCS superconductor is ~ TVAq , 
where ./V is the density of states at the Fermi level. Since 
the coherence length is vf / Ao one finds that the core en- 
ergy of a vortex is ~ e/r, the Fermi energy, which is much 
larger than the BCS superconducting T c . This is why 
the superconducting transition of a weakly coupled su- 
perconductor is essentially mean-field in character, with 
an unobservably narrow critical region. It seems then 
that a competing order developing in the vortex core is 
necessary to bring the core energy down, so that vortices 
may proliferate at a low T c . That this, however, does 
not necessarily follow, one may see by realizing that the 
above core energy (~ ep) is nothing but the bare stiff- 
ness for the phase fluctuations. The role of the compet- 
ing order parameter would then be to reduce this stiffness 
with underdoping. But as we argued, this seems perfectly 
possible without any ordering, simply from the repulsion 
which suppresses phase coherence even at T = 0. As 
long as the T = stiffness is small, the BKT transition 
temperature will also be small |2^. The problem with 
the argument in favor of the competing order is that it 
ignores possible quantum fluctuations arising from inter- 
actions, which are on the other hand, the main point of 
the present work. 
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VII. APPENDIX: DERIVATION OF THE 
BOSONIC MEAN-FIELD ENERGY 

Here we present the details of the derivation of the Eq. 
(23) for the bosonic mean-field energy. Within the local 
approximation discussed right above Eq. (22), the parti- 
tion function for the effective bosonic system becomes 

Z B = J D[^*,y}Dp-e~fo dTLB j (35) 

with the bosonic Lagrangian 

L B {T)=Y,^*{xiT)d T ^{x,T)- (36) 

X 

t B J2 *>,t)*(x',t) 

(x,x f ) 

+ U'J2\*(^T)\ 2 (\y(x,T)\ 2 -l) 

X 

+C/^|*(x,r)| 2 p-(^,r) 

X 

t)\ 2 + n-p-(x, t)) 

X 

+U'J2p-&t)(p-(x,t)-1), 

x 

where we have written the hard core interaction U' 
(U' — ► oo ) explicitly. This way one can still write the 
standard coherent state representation of the partition 
function. One can then derive Eq. (23) as follows. First, 
decouple the hopping term by introducing the Hubbard- 
Stratonovich complex field <!>(x,t): 

J ^ E<.,./) **(«.r)*(-'.r) K J D[ <p %m) 

-J <M^ x x , **(^,T)(*B 1 )x.x'*( a: ! ' r )+X]x *( 2; ^' r )**( 2; ^ T )+ C - C -] 

Next, calculate the partition function in the saddle point 
approximation. Assuming that at the saddle point 
$(x,t) = and p-(x,r) = p~[x), i. e. the static 

configuration, in the hard-core limit U' — > oo the parti- 
tion function becomes 

Zb mf = e -*E^, ••W(*i 1 WC) Tre -5 ¥t> (38) 

where the mean-field Hamiltonian is 

Hmf - £(*(aO&t(aO + **(x)b{x)) ^39) 

X 

5> + &t(x)6(x) + fi-p-(x)) +Uj2P(x)K*)P-(x), 

X X 



9 



and the trace is to be evaluated over the four states as can be easily computed. Finally, rescaling <f> — > $/2 leads 
described below Eq. (23). Since Hmf is local it factorizes to the mean-field free energy as displayed in Eq. (23). 
into a product over the lattice sites, and thus the trace 
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